Magnetic response enhancement via electrically induced magnetic moments 
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The realization of negative refraction in atomic gases requires a strong magnetic response of 
the atoms. Current proposals for such systems achieve an enhancement of the magnetic response 
by a suitable laser field configuration, but still rely on high gas densities. Thus further progress 
is desirable, and this requires an understanding of the precise mechanism for the enhancement. 
Therefore, here we study the magnetic and electric response to a probe field interacting with three- 
level atoms in ladder configuration. In our first model, the three transitions are driven by a control 
field and the electric and magnetic component of the probe field, giving rise to a closed interaction 
loop. In a reference model, the coherent driving is replaced by an incoherent pump field. A time- 
dependent analysis of the closed-loop system enables us to identify the different contributions to the 
medium response. A comparison with the reference system then allows one to identify the physical 
mechanism that leads to the enhancement. It is found that the enhancement occurs at so-called 
multiphoton resonance by a scattering of the coupling field and the electric probe field mode into the 
magnetic probe field mode. Based on these results, conditions for the enhancement are discussed. 

PACS numbers: 42.50.Gy, 42.65.Sf, 42.65.An, 32.80.Wr 



I. INTRODUCTION 



As Veselago pointed out in 1968, materials having both 
negative permittivity and permeability can acquire a neg- 
ative index of refraction [l|, Q • Such negative index ma- 
terials are also called "left-handed" , since the electric 
(E) and magnetic (H) components of an electromagnetic 
wave travelling trough a negatively refracting medium 
and its wave vector form a left-handed coordinate system. 
These materials offer promising applications 0, [H, 0|, 
such as the possibility to overcome the diffraction limit 
with a negatively refracting, perfect lens, as proposed 
by Pendry [f|. Therefore it is not surprising that in the 
recent past, left-handed materials and negative refrac- 
tion have been studied intensely. These efforts have been 
fueled by a multitude of successful experimental demon- 
strations of negative refraction and related effects, mostly 
relying on metamaterials d, 0, H, H, EH H3 ■ Metamate- 
rials are artificial structures with feature size below the 
incident radiation wavelength that allow to control the 
electromagnetic response to a great extent and at the 
same time appear as a bulk medium to the incident ra- 
diation. 

A different ansatz is the quantum optical realization of 
negative refraction in dense atomic gases [HI, [l3|, 0, [ll| . 
Negative refraction, however, typically requires both elec- 
tric and magnetic response at the probe field frequency, 
which is hampered by the fact that usually the coupling of 
the magnetic component of a probe laser field to a mag- 
netic dipole transition in atoms is strongly suppressed. 
A simple order-of-magnitude estimate shows that the 
suppression factor is proportional to two powers of the 



'Electronic address: 
t Electronic address: 



bastian.jungnitsch@mpi-hd.mpg.de 
joerg.evers@mpi-hd.mpg.de 



finestructure constant, a 2 ~ 137~ 2 . Thus, additional 
effort is required in order to enhance the magnetic re- 
sponse. In the literature, a number of schemes that allow 
to achieve a high positive index of refraction with small 
absorption have been proposed. They are based on a suit- 
able modification of the electric response of the medium 
(see, e.g., (HI). The enhancement, however, is typically 
too small such that a direct transfer of these ideas to 
magnetic transitions is not straightforward. Thus, the 
schemes suggested so far for negative refraction rely on 
a different mechanism that is related to an enhancement 
via chirality [H, EH, EE E3l • The medium is such that the 
magnetic response is influenced by both the electric and 
the magnetic component of the probe field and analo- 
gously for the electric response. A first interpretation for 
the case of atomic systems has been given in [l2j , which 
however did not consider the coupling of the magnetic 
probe field component to the atomic system in deriving 
the induced magnetic dipole moment, and which focusses 
on a certain resonance case for the applied fields and the 
employed level structure. Also, the system considered 
there only allows for enhancement at a single frequency, 
while an enhancement at a range of probe field frequen- 
cies was reported in subsequent work [lit [l5| . Thus 
better insight is desirable, not least since it might lead 
to further enhancement of the magnetic response which 
would significantly simplify the theoretical and experi- 
mental study of negative refraction in atomic gases. 

Motivated by this, here we study in detail the en- 
hancement mechanism that is at the heart of current 
schemes to achieve negative refraction in atomic gases. 
For this, we revisit the three-level ladder-type system 
studied in [l2j , where one transition is driven by a coher- 
ent control field and the other two transitions couple to 
the magnetic and electric component of a probe field. In 
contrast to previous studies, we apply a time-dependent 
analysis of the medium response that enables us to di- 



rectly identify the various physical processes contribut- 
ing to the medium response. These results are compared 
to a reference system that is obtained by replacing the 
coherent driving field by an incoherent pumping field. 
We find that the enhancement of the magnetic response 
occurs since the used level schemes are so-called closed- 
loop media. In these systems, the laser fields applied are 
such that they form a closed interaction loop. We iden- 
tify a scattering of the coupling field and of the electric 
probe field into the magnetic probe field component as 
the mechanism responsible for the enhancement of the 
response and provide conditions for this process to take 
place. It is found that the so-called multiphoton reso- 
nance condition must be fulfilled for the enhancement to 
be present. In the studied three-level system, the con- 
dition is satisfied only at a single probe field frequency. 
But in larger level schemes, the laser fields can be applied 
in such a way that the enhancement works at arbitrary 
frequencies of the probe field. 

The paper is organized as follows: In Sec. [Til we present 
our model systems and derive the equations of motion as 
well as expressions for the medium response coefficients. 
In Sec. IIII1 we solve the equations of motion, both in the 
time-dependent case for the closed-loop configuration, for 
the time-independent case at multiphoton resonance and 
in the incoherently pumped system. Using these results, 
in Sec. IIVI we compare the different systems and finally 
identify the enhancement mechanism. Sec. [V] discusses 
and summarizes the results. 

II. THEORETICAL CONSIDERATIONS 
A. Model 

We start by writing down the applied electromagnetic 
fields as shown in Fig. [TJ Since we treat both systems 
semi-classically, we have 

E(f, t) = Epirje^e-^"* + EJtfe** e^*"* 

+ c.c. , (la) 

B{r, t) = Sp^e^e-**'* + B c (r)e^ e"^* 

+ c.c. , (lb) 

where the subindex p [c] refers to the probe [control] field. 
Further, E p (r) = £ p e p e' lkpr , where £ p is the electric field 
amplitude, e p the unit polarization vector of the electric 
component of the probe field, k p is the probe field's wave 
vector, Up its frequency and (j) its absolute phase. Anal- 
ogously, E c (r) — £ c e c e lkaV , where £ c is the control field 
amplitude and e c its unit polarization vector. k c is the 
wave vector of the control field, oj c its frequency and ip 
the control field's total phase. The magnetic probe field 
component is defined analogously as B p (r) — B p b p e lkpr . 
Note that the magnetic probe field component unit po- 
larization vector is b p — k x e p with unit propagation 
direction vector k = k p /k p . 
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FIG. 1: (Color online) (a) The three-level system driven by 
coherent fields in loop configuration. The probe field compo- 
nents are denoted by red solid double arrows, the coupling 
field by purple solid double arrows. Spontaneous emission is 
indicated by the wiggly green arrows. The transition |1) <-> |2) 
couples to the magnetic component, while transition 2) <-> |3) 
couples to the electric component of the same probe field, (b) 
Reference system obtained by replacing the coherent control 
field by an incoherent, bi-directional pumping, indicated by 
the green dashed double arrow. 

In rotating-wave and dipole approximation, we arrive 
at the Hamiltonian 

H =H + H I , (2a) 

3 

i=i 

Hi = -h (0 2 ie-^ t |2)(l| + ft 32 e-^*|3)(2| 

+ r2 3 ie-^ t |3)(l|+H.c.) . (2c) 

Note that in the rotating wave approximation, the mag- 
netic control field component can be neglected, because it 
does not couple near-resonantly to a magnetic transition. 
In Eqs. the energy of state \i) is hui (i € {1,2,3}) 
and the Rabi frequencies are defined as 

n 21 =e^B p (r)fI 21 /h, (3a) 
n 32 =e i '»E p (r)d 32 /h, (3b) 
fi 3 i =e^E c {f)d 31 /h. (3c) 

The electric dipole moments are defined as d 32 — (3|d|2) 
and G?3i = (3|rf|l). Analogously, the magnetic dipole mo- 
ment is fi2i = (2|/2|1), with fli2 = Here, d and jl are 
the electric and magnetic dipole operator, respectively, 
and dij = d* t , = and fiy = 

B. Equations of motion 

In this section, we derive the equations of motion for 
a general system that contains both systems of inter- 
est shown in Fig. [1] as special cases. To this end, 
we consider a three-level system that - on transition 
1) «-> 1 2) - combines coherent pumping by a control 
field and bi-directional incoherent pumping. Since this 
system is a closed-loop system, in general there is no 
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stationary state in the long-time limit and the Hamil- 
tonian necessaril y h as an intrinsic explicit time depen- 
dence [H, d MM, M, 0, S i3, H H m ■ 

To simplify this time dependence, we apply the transfor- 
mation 



(4) 



where X = Ai|l)(l| + A 2 |2)(2|. Here, we chose the no- 
tations Ai = — lj\ — oj c and A 2 = 0J3 — u>2 — u p for 
the detunings. By applying this transformation, we ar- 
rive at the following equations of motions, if we include 
spontaneous decay in the Born-Markov approximation: 



d 

d „ 
dl 912 



dl Ql3 



ngii + 73022 + (71 + n) 033 + ie 

A '^21012 + ^13031 - ^31013 , 



-iAt, 



ie 



Sll2021 

(5a) 

012 - «^32013 
~ lA *^i2 (011 ~ £22) + iftl3f?32 , (5b) 



1 



i(A-A 3 ) + -(n+7 3 ) 



i(A 
ifii 3 (gu 



A 2 - A 3 ) + i (2n + 71 + 72) 



013 



dt 



022 



Q33) - 1VL23Q12 + ie 
73022 + 72033 - *e~* A *Sli2f?2i 

iAt 



ie 



£^21012 ~ 2^32023 + 2^23032 



iAt ~ 

"12 023 j 

(5c) 
(5d) 



d . 



-jA 2 

2^13021 



2 v 1 



+ 7i + 72 + 73) 



023 



ie 



•AI 



^21013 + i^23 (033 - Q22) 



033 =1 ~ 011 - 022 



(5e) 
(5f) 



Here, Qij £ {1,2,3}) is the density matrix in the 
interaction picture obtained by transformation of ac- 
cording to Eq. (0| . One can see that this transfor- 
mation simplifies the explicit time dependence on the 
right hand side of the equations of motion to factors 
of e ±lAt in front of the weak magnetic probe field Rabi 
frequency f2 2 i or O12, respectively. 7, are spontaneous 
emission rates on the different transitions. Also, we intro- 
duced the detuning on the magnetic probe field transition 
A3 = ll>2 — u\ — Up, as well as the so-called multiphoton 
detuning 



A = A 2 + A 3 - Ai , 
which for the current system evaluates to 
A = u)n — 2u r . 



(6) 



(7) 



By setting the incoherent pumping rate r± — in 
Eqs. (O, we arrive at the equations of motion for the 
system in Fig. [TJa). Setting £! 31 = £! 13 = A 3 = yields 
the equations of motion for the incoherently pumped sys- 
tem shown in Fig. HJb). 



C. Electric and magnetic responses 

Since our aim is to study the magnetic and electric re- 
sponses, we require an expression for them in terms of 
the density matrix elements governed by Eqs. ©. We 
will find such a relation in this subsection. For this, it is 
important to note that electric fields can not only induce 
electric polarization, but also magnetization [3, 53] • 
Similarly, magnetic fields can induce both magnetization 
and polarization. 

For definitiveness, in the following we specialize to 
a circularly polarized (<r + ) probe field and probe field 
propagation in z direction, since one then obtains for 
the probe field polarization vectors b p = —ie p , i.e., they 
are parallel. Therefore, the tensorial structure of the re- 
sponse coefficients in the macroscopic polarization P and 
magnetization M simplifies considerably 



P(f,t) 



M(r, t) 



£eh(t) H(f,t-r) dr 



+ £0 / Xe(r) E(r,t - t) dr , (8a) 

J — OO 

/ Zhe(t) E{r,t-r)dT 

c Mo J-00 

poo 

+ Xm (r)H(r,t-T)dT, (8b) 



and the electric and magnetic susceptibility Xe and 
Xm and the chirality coefficients £,he and £eh become 
scalars. While the susceptibilities determine the electric 
[magnetic] response to the electric [magnetic] probe field 
component, the chiralities or cross-terms determine the 
magnetic response to the electric probe field component 
and vice versa. Here, c is the vacuum speed of light, and 
Ho and So are the vacuum permeability and permittivity. 
Note that, with the notation of Eqs. {§]), the refractive 
index in Fourier space is given as [14| 



»U') = \/ i(ui)jj,(u>) - - £eh(~<-') + O/i: U\) 



1 

+ 2 



£,eh{u) -£he{w) 



where 



eM = XeM + 1 ■ 
MM = XmM + 1 ■ 



(9) 



(10a) 
(10b) 



are the permittivity s(uj) and the permeability jx{ui) in 

Fourier space. Also, £,eh{u), £he(w), Zm(w) and x e (cj) 
are the Fourier transformed response coefficients of the 
corresponding time domain quantities in Eqs. l[8"]). They 
satisfy Xe{^)* = Xe{— and similar for the other sus- 
ceptibilities and the chiralities. For vanishing chirality 
coefficients, Eq. © reduces to the well-known relation 
between the refractive index and the permittivity and 
permeability for non-chiral media (32j. 
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Let us now continue to find an expression of the re- 
sponse coefficients in terms of the density matrix ele- 
ments. Plugging Eqs. (TT]) into Eqs. (0), we arrive at 

p = j E H(u P ) g (r y(0-^) 

C/i 

c^o 

+ e Xe{u c )E c {r)eW-"^ + c.c. , (11a) 

M = ^K) ^g (f)e i(0-o> pt ) 
Hoc 



-X m (uj p )B p (r)e 1 ^-^ 



Mo 

+ — Xrni^M^e^-^ + C.C. . (lib) 

Mo 

Here, we have used B = fioH, which holds since B is 
an external field. Note that the Rabi frequencies in the 
equations of motion contain local fields, whereas Eqs. ([8]) 
contain external fields However, local field effects 

are not considered here and therefore we do not distin- 
guish between external and local fields, as is valid for 
moderate particle densities. 

The polarization is also given by P = Np, and the 
magnetization by M — Nm, where N is the particle 
density and p and m are the mean polarization and mag- 
netization per atom [33| . We can express the mean po- 
larization as p — Tr(gd) and the mean magnetization as 
to = Tr(gfl). These traces can be written in terms of the 
transformed density matrix elements s {1, 2, 3}) 

which are given as solutions of the equations of motion 
([5]) . If the light travels through the medium over a macro- 
scopic distance, then the magnetic and electric response 
is determined by the part of the medium response that os- 
cillates in phase with the probe field. In order to identify 
the relevant parts, we apply another unitary transforma- 
tion into a system oscillating in phase with the probe 
field. We denote the density matrix in this frame by 

QlJ e {1,2,3}). 

In the incoherently pumped system of Fig. QJb), the 
Hamiltonian is time-independent, and we find 

Q21 = 021 , (12a) 

§32 = 032 ■ (12b) 

On the other hand, in the closed- loop system of Fig.QJa), 
the coherences £32 and £21 are given by 



021 
932 



„iu> v t 

e p Q21 



e p ' Q21 



e p Q32 = Q32 ■ 



(13a) 
(13b) 



One can see that at multiphoton resonance A = <^> 
uj c = ujp/2 [see Eq. (J7J], the two reference frames de- 
noted by pij and coincide for the two coherences in 



Eq. p3)l . similar to the case of incoherent pumping. We 
will investigate the multiphoton resonance case further 
in Sec. IIIIA2I 

We now proceed with the evaluation of the response 
coefficients. Keeping only terms oscillating in phase with 
the probe field, we arrive at 



Xe 



N 

s h 



^32 032 ^ ! 



N ^0 2 ^(1,0) 
fc M21 021 



~ . Nc^o , 40,1) 

£he = -1 — z — "32 M21 e 02i 



EH 



h 



(14a) 
(14b) 
(14c) 



d 3 2 M21 e- 4 *^ , (14d) 



where the g\j ' and g\j : are expansion coefficients in a 
Taylor series of Qij (i, j S {1, 2, 3}) in terms of O32 and 

032 — 032 + 032 "32 + 032 J '21 

+ o(n^ I n§ 2 ,n 2 in 32 ), (i5a) 

021 — 021 + 021 "32 + £> 2 1 "21 

+ 0(^,0^,021^32). (15b) 

In Eqs. (fT5|) . we call ggg^Oaa and £? 2 i'°^2i the "di- 
rect terms" , since they correspond to the susceptibilities, 

while £> 3 2'°''02i and g^^fl^ are denoted "cross terms" 
as they give rise to the chiralities. Also, in Eqs. (Tj"4")) . we 
have introduced the relative phase, 

$ = 032 - ^21 , (16) 
between the scalar dipole moments which we write as 

d 32 e p = d 32 , (17) 



M21&P — M21 1 



(18) 



Eqs. (|T4|) are the desired relation between the density 
matrix elements and the coefficients that determine the 
magnetic response. These can now be used in Eqs. (1111) in 
order to determine the contribution of the various pro- 
cesses to the polarization and magnetization. Keeping 
only the terms relevant to the probe field response in 
phase with the probe field frequency, we find 



N 



-i* -(1,0) 



+ d 32 £p P32 ) + cc. , 



N 



2 R 

M21 °P 021 



C.C. 



(19a) 



(19b) 



III. ANALYTICAL RESULTS 

We next solve the time-dependent equations of mo- 
tion (JSJ) for both of our systems. First, we consider the 
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closed-loop system for arbitrary multiphoton detuning A 
and derive an expression for the coherences to first order 
in the magnetic and electric probe field Rabi frequencies 
0,21 and f^32. Then, we consider the special cases A/0 
and A = 0. Finally, we solve the incoherently pumped 
system. 

A. Closed-loop system 

The equations of motion for the closed-loop system are 
obtained from Eqs. {!]), if we set n — 0. We define the 
vector R containing all density matrix elements, 

R = (qu,Qu, §13, 021, 022, 023, 031, 032) T ■ (20) 

The equations of motion (JSJ) can be rewritten in terms of 
R as 

—R = MR + T, . (21) 
at 

Here, we have eliminated 033 via the trace condition 
Tr(g) = 1, which is the reason for the appearance of 
the constant term S in Eq. (JUJ). We proceed by split- 
ting both M and R up into terms with different time- 
dependencies as follows: 

M = M + MiQ 2 ie tAt + Af_i«i 2 e~ iAt , (22) 
E = E + Si« 2 ie lAt + £_itt 12 e- iAt , (23) 

where Mk and (k € { — 1, 0, 1}) are time-independent. 

According to Floquet's theorem [34[, the solution of R 
has only contributions oscillating with frequencies that 
are integer multiples of A. Since terms oscillating at 
higher frequencies are suppressed by powers of the mag- 
netic probe field Rabi frequency |^2i|, we expand R to 
first order in this Rabi frequency, and work with the 
ansatz 

R = R + Rin 2 ie lAt + R^iQi 2 e- iAt + 0(\Q 12 \ 2 ) . (24) 



From Eqs. (|2HI24[) . a comparison of coefficients yields 

Ro = -Mq-^q , (25a) 
Ri = -(M - iA)-\MiR + Ei) , (25b) 
R-i = -(M + iA) _1 (M_i^o + E-i) • (25c) 

Since the density matrix element 021 oscillating in phase 
with the probe field and 021 are related as 

021- e-^-^hi 

= e-^- 2 "*)* [R U + [R1USI21 

+ e-WK-awpJt [fl_i] 4 fii 2 t (26) 

one can determine [Ri]^ as the part of 021 oscillating in 
phase with the probe field, where the index 4 denotes the 
fourth component of Ri. Likewise, the relevant part of 
032 can be identified with [i?o]s- 

Note that Ro, J?_i and Ri are independent of f^i, 
but do depend on the electric probe field Rabi frequency 
^32, which we did not take into account so far. Since 
we are only interested in the linear magnetic and electric 
response, we still have to expand the appropriate com- 
ponents of R in fi 32 . We obtain 

21 =[i? 1 ] 4 fi 2 i + [i?o]4e- l ^- 2 ^ t 

+ o(fi^,n| 2 ,n 21 n 3 2), (27a) 

032 =[R-i]s e-^'-^nia + [Roh 

+ o(n^,o 2 2! n 21 n 32 ) ! (27b) 

where [i?ok <x f2 2 3 and [Rq]s oc f2 3 2, and higher orders 
of fl 32 have been neglected. With Eqs. ([zTS]l and P7|). an 
explicit evaluation yields 



021 =77 



032 



8A? 73 + 4iAf 73 (2iA 3 + 7s ) + A 3 [8|fi 3 i| 2 (72 



73 j 



2r 



4j0 3 i | 2 + (2A 3 - 173) (2Ai - 2A 3 + 17.) 

2Ai [— 4|Q 3 i| 2 (72 ~ 2 73 ) +T\-i [(4|0 3 i | 2 72 - T) 7. - 4|» 3 i| 2 7 3 2 ] 
4|fi 3 i | 2 + (2A 3 - z 73 ) (2Ai - 2A 3 + i ls ) 



4 4|ft 31 | 2 (72 - 73) + (2Ai + i( 7 i + 72)) 73 (-2A 2 - i 7s ) 



B 



4|r> 3 i| 2 + [2i (Ax - A 2 ) + 73] (-2«A 2 + 7s ) 



-i(uc-2u> P )t Q 31 Q 



23 , 



4 4|0 3 i| 2 (-72 + 73) + [2Ai +i (71 + 72)] 73 (2Ai - 2A 3 - i 7s ) 



D 



4|f» 3 i| 2 + (2A 3 + z 73 ) (2Ai - 2A 3 - i ls ) 

,8 [2(-Ai + A 2 )72 + H3(2iA 2 +7i+73)] 
fl 4|0 3i | 2 + (-2iAi + 2 2 A 2 + 73) (2zA 2 + 7s ) 



n 



32 



(28a) 



(28b) 
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where 

B = 4A^73 + T + 4|0 31 1 2 (72 + 273) , 

7s = 7i + 72 + 73 > 

T = (71 + J2) 2 73 • 

Note that in Eq. (|28|) , the multiphoton detuning A has 
been replaced using Eq. ([6|). Also, the contributions 
without the explicit time dependence via an exponential 
factor are the direct terms, while the other parts are cross 
terms. 



1. Non-zero multiphoton detuning (A 7^ 0) 

In the case of non-zero multiphoton detuning, A =/= 
<^> 2u) p ^ w c , only the direct terms in Eqs. (|28p oscil- 
late in phase with the probe field. Therefore, only these 
terms contribute to the coherences and thus to the mag- 
netic and electric response. Then, the cross terms in 
Eqs. (fT4")l vanish and therefore the chirality coefficients 



vanish, £eh — £,he = 0. This means that the polariza- 
tion [magnetization] is entirely determined by the elec- 
tric [magnetic] probe field component. It will turn out in 
Sec. lIII Bl that this case is comparable to the incoherently 
pumped system, in which there arc no cross terms cither. 



2. Multiphoton resonance (A = 0) 

We now focus on the case of multiphoton resonance, 
i.e. A = or u p = ui c /2. Hence, Eqs. ([5]) become 
time-independent and we can now solve the linear sys- 
tem Eq. (f2"Tj) for a time-independent steady-state solu- 
tion of R using Jj.R = 0. Now, all terms in Eq. (j2"4"|) 
contribute, apart from terms that contain SI32 in higher 
than first order. From Eqs. (fT3|) we also find that the 
simplified relations between the two considered reference 
frames Eqs. (fT2|) hold as in the case of incoherent driving. 

We again neglect terms of higher order in the probe 
field Rabi frequency and arrive at 



2i 

hi = - 77-77 {SiA:^ + 4A 2 (4zA 3 - 7s ) 73 - ^73 + r) 7s + 4|ft 3 i | 2 [72 (71 + 72) 



C+D 

(2iA 3 + 72 ) 73 - 7 f] + 2iA 2 [-4|ft 3 i | 2 (72 - 2 73 ) + T + 4 73 A 3 (A 3 + 17,)] } n 21 

{4|f>3i| 2 (72 - 73)-[2(A2 + A 3 )-H(7i+72)]73(2A2+i7 s )} ^23, (30a) 



4 ^31 f , „ J 

C+D 
4fl 



032 = { 4 l^3i| 2 (72 - 73) - [2 (A 2 + A 3 ) + i (71 + 72)] 73 (2A 2 - i ls )} Oi 



C-D 



■fi 3 i I 2 [2i A372 + 73 (2iA 2 + 71 + 7s)] «32 , (30b) 



where B. Incoherently pumped system 



C± =4|r» 3 i| 2 + (±2iA 3 + 73)(T2 l A2+7.) , 
D = 4 (A 2 + A3) 2 73 + T + 4|f> 3 i| 2 (72 + 2 73 ) , 

r = (71 + 72) 2 73 1 

Is = 7l + 72 + 73 ■ 

The control field detuning Ai has been eliminated using 
the relation Ai = A2 + A3 which follows from Eq. ([6]) in 
the case of A = 0. 

We see that both the direct terms and the cross terms 
contribute to the coherences and thus to the magnetic 
and electric response in this case. Therefore, the chirality 
coefficients are non-zero in the resonance case A = 0, see 
Eqs. (O. 



We now consider the electric and magnetic response 
in the incoherently pumped system shown in Fig. [TJb). 
It will serve us as a reference in a comparison to the 
results of the closed-loop system in order to determine the 
exact origin of the various contributions to the medium 
response. 

In this system, the transformed Hamiltonian is time- 
independent. Hence, we can solve for the steady-state 
solution just as in the case of multiphoton resonance in 
Sec. lIII A~2] The equations of motion follow from Eqs. §5§ 
with Ai = and SI31 = ^13 = 0. Instead of the coherent 
coupling field, they include an incoherent, bi-directional 
pump rate r\ on transition |I) <-> |3). 

Using Eqs. (|I2p. up to first order in the probe field we 
find for the coherences in a reference frame oscillating in 
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phase with the probe field: 

2i[n(7 3 - 72) + (71 + 72)73] 



02 1 



032 



n 2 



(2iA 3 + n + 73)[(7i + 72)73 + ^1(72 + 273)] 
+ 0(^,^,^032), (32a) 

2ir x (7 2 - 73) 



32 



(2iA 2 + 7s + ?*i ) [(71 + 72)73 + n(72 + 273)] 
+ 0(^,0^,0^032). (32b) 

In this case, the electric [magnetic] probe transition co- 
herence is determined by the electric [magnetic] probe 
field component, and no cross-terms appear. 
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IV. COMPARISON OF THE TWO SYSTEMS 

We will now proceed with a comparison of the two 
systems. To this end, we will discuss the expansion co- 
efficients £32'^ and £21 i ( a 'k e {0, 1}) in Eqs. (fl~5|) . 
since they determine the magnetic and electric response 
according to Eqs. (fl9|) . 

It will turn out that the direct terms are similar in 
many regards in both systems. These terms describe the 
establishment of polarization [magnetization] due to the 
electric [magnetic] probe field component. But crucial 
differences are found for the cross terms, which only ap- 
pear in the closed-loop system for zero multiphoton de- 
tuning. These terms characterize the polarization [mag- 
netization] due to the magnetic [electric] probe field com- 
ponent. 



The direct terms 



From Sees. IIIIAI1 MI A 21 and [ITTBl it is clear that di- 
rect terms appear both in the loop system and in the 
system with incoherent pumping. In Fig. [51 the cor- 
responding expansion coefficients for the two cases are 
compared. Since transitions |1) <-> |2) and |1) <-> |3) are 
electrically dipole-allowed, whereas transition |2) |3) 
is magnetically dipole-allowed, we choose the decay rates 
as 71 =7,72 = 7,73 = a 2 7- For the incoherent case, 
the pump rate is set to 7*1 = 7, whereas in the closed- 
loop configuration, the coherent pump field is set to 
Ai = 0, — 7. All expansion coefficients are plot- 
ted against the respective probe field detunings A2 or 
A 3 . 

We find that in many respects, the direct terms of both 
systems behave similarly. Apart from the AC Stark split- 
ting in the closed-loop system, the general structure is 
comparable and in particular the magnitude of the coef- 
ficients is of similar order in both systems. 

The similarities become more apparent when consid- 
ering the dependence of the direct terms on the coher- 
ent pump rate f^i and the incoherent pump rate n, as 
shown in Fig. [3J For this figure, the direct terms of the 
loop system are evaluated at a particular detuning A2 



FIG. 2: (Color online) Real (red solid curve) and imaginary 
part (blue dashed curve) of the direct terms in the magnetic 
and electric susceptibility, respectively. The top row shows 
g2x 1 the bottom row £32 • (a) Closed-loop system and (b) 
incoherently pumped system as shown in Fig. Q] 



or A3, respectively, at which the absolute value of the 
imaginary part becomes maximal. At the same point, 
the real part vanishes. For the incoherently pumped sys- 
tem, the corresponding maxima always occur at A3 = 
or A2 = 0, respectively. This approach allows to compare 
the two systems independently of the AC Stark splitting 
appearing in the closed loop system only. Due to the 
growing splitting with increasing IO31 1, a comparison at 
a fixed detuning would not be meaningful. It can be seen 
from Fig. [3] that both systems show a qualitatively sim- 
ilar dependence on the pumping strength. The shown 
imaginary part of the coherences characterizes the ab- 
sorptive behavior of our systems: positive values stand 
for absorption and negative values for amplification of 
the magnetic or electric probe field component. 

Interestingly, in both systems it is possible to choose 
the respective pump rate in such a way that 021 vanishes 
at all frequencies. This is the case at the roots of the 
dashed lines in Fig. [3] It turns out that at these points, 
the populations of states |1) and |2) are the same such 
that the magnetic probe field component can traverse the 
medium without attenuation and without experiencing 
diffraction. 

For the interpretation of Fig. [3J we calculate the co- 
herences in terms of the (zeroth order) populations for 
arbitrary detuning. In the case of incoherent pumping, 
we obtain 



021 



032 



2^21 



-(0) -(0) 

2A 3 -i(n+7 3 ) ' 

-(0) _ -(0) 
022 033 



' 2A 2 - i(n + 71 + 72 + 73) 
The zeroth order populations are 



.(0) (n + 71 + 72)73 
011 = 



C 



(33a) 
(33b) 

(34a) 
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FIG. 3: (Color online) Dependence of the expansion coeffi- 
cients §32 in the closed loop (black solid line) and in the 
incoherently pumped system (green long dashed line) as well 
as p 2 i'°' f° r the closed-loop (blue dash-dotted) and the in- 
coherently pumped system (red dashed) on the strength of 
the control field and the incoherent pump rate, respectively. 
Shown are only the imaginary parts at the maximum of the 
absolute value of the coefficients. For the incoherent config- 
uration, the maximum of f?2i'°' is always at A3 = 0, and the 
maximum of §32 is at A2 = 0. The coefficients are scaled by 
7 to obtain a unitless quantity and §21'°' is scaled by a factor 
of 10" 4 . In this fi gure, negative values indicate amplification, 
positive values absorption. 



40) 
022 = 

40) 

033 = 



7-172 

c 

7-173 

C 



(34b) 
(34c) 



show that 



^ "(0) j 

< Q22 and 



where C = 7772 + 27773 + 7173 + 7273. 

Note that in our system, 73 <C 72 due to the differ- 
ent multipolarity of the transitions. Therefore, for small 
pump rates r\, from Eqs. (f3"4"| one finds g^ > g^2 > 633 ■ 
As a result, both probe transitions are absorptive - al- 
though for P21 only in a very small range of 77 as com- 
pared to 7. 

For large 77, Eqs 

033 < 822- I n f &c t> Oil ~ Q33 f° r r i ^ 7- Hence, the 
magnetic probe transition becomes amplifying, whereas 
the electric transition maintains its absorptive character. 
As expected, from Eqs. (|33p . it follows directly that a 
population inversion causes amplification. 

As a side note, we would like to mention the fact that 
Eqs. (f3"3"| imply vanishing £32 for 72 =73. This is due 
to the fact that from Eqs. (j3~4"|) , one then finds g§} = 
ggg* . However, this case is not of relevance for the current 
analysis, since 73 <C 72- 

Let us now examine the behavior of the closed-loop 
system with regard to a change of ^31 1. In this case, the 
coherences are given by 



(4|0 31 | 2 + ^2)^3 



(0) 



021 



2Qo 



932 2 n ' 2 (fW + i^ 



while the populations obey 



-(0) 

0n 



40) 
Q22 



»(o) 

033 



(35b) 

(36a) 
(36b) 
(36c) 

r i72|"3ir 1" °73|"3i|" 
G {1, ... 7}) are coefficients independent 



73(^6 + 4|ft 3 i| 2 ) 


K 7 - 


f 472^31 1 2 + 873 


^3l| 2 




472 1^31 1 2 




K 7 - 


f 4 7 2 10 3 1 1 2 + 873 


^3l| 2 




4731^31 1 2 




K 7 - 


f 472^31 1 2 + 873 


^3l| 2 



where the Ki , 
of f2 3 i . 

Again, the coherence Qij is determined by the dif- 
ference of the populations of grp and g^ , € 
{(2, 1), (3, 2)]. For small [O31 1 , the term proportional to 

033 ~ 011 can be neglected. For large |f2 3 i | , this term out- 
weighs the others at first sight. However, for |£l3i| — * 00 
and arbitrary, but fixed detunings, one finds 



d (o) 
0ii 



40) 

022 



03^ 



1 



21 + 2 

73 



22a + 1 ' 

72 



21 + 2 

73 



(37a) 
(37b) 
(37c) 



Therefore, g{i 



such that also in this case, g*ij 

"(0) 

0jj 



(35a) 



-(0) 

033 

is determined by gf^ 

Eqs. (|3~7| also imply that for a strong control field, 
both probe transitions show opposite absorptive behavior 
(see Fig. [3]): absorption in case of the upper probe tran- 
sition, and amplification for the lower one. The reason 
is that due to the small decay rate 73, most population 
is trapped in state |2), such that the relevant population 
differences in Eqs. (|3"5"|) have opposite sign. For small 

|^3i |, Eqs. (|3"6"|) yield g$ > g§} > £33 , which explains 
the absorptive properties in Fig. [3] However, the first 
inequality is not obvious and can only be deduced by an 
exact knowledge of K 6 . 

Before we come to the discussion of the cross terms, we 
would like to note that, of course, there are also differ- 
ences between the direct terms of the two systems. We 
already discussed the AC Stark shift that occurs only in 
the loop system. But the closed-loop system also offers 
more degrees of freedom than the incoherently pumped 
system, most importantly Ai. The role of the control 
field detuning Ai is shown in Fig. 2] In this figure, we 
plot the expansion coefficients of the direct terms in the 
coherence over the probe field detunings A 2 and A 3 for 
Ai = 27. The detuning Ai essentially determines the 
position of one of the maxima of the imaginary part of 
the response function. Increasing Ai moves the peak to 
higher frequencies, while decreasing it moves it to lower 
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A 2 / 7 

FIG. 4: (Color online) Real (red solid line) and imaginary part 
(blue dashed curve) of the expansion coefficients (a) §21 an d 
(b) P32' 1 ' in Eqs. (|15[1 . The curves are drawn for Ai = 27 in 
the closed-loop system. The coefficients are scaled by 7. The 
parameters are as in Fig. [21 a), except for the non- vanishing 
detuning of the control field. 



frequencies. For example, in Fig. 31 the corresponding 
maxima can be seen close to A2 = 27. 



B. The cross terms 

Let us now come to the most important difference be- 
tween the two systems: The coherences of the closed-loop 
system have cross terms, while the coherences of the in- 
coherently pumped system do not. However, as found in 
Sec. IIII A) the cross terms only contribute to the mag- 
netic and electric response for A = 0, i.e. for lu p = uj c /2. 
The following comparison serves as basis for the conclu- 
sion which will be drawn in the discussion section regard- 
ing the enhancement of the magnetic response. 

The mathematical origin of the cross terms can be 
identified from the derivation of the coherences in Sec- 
tion IIII Al The relevant response of the system to the 
probe field is given by the contributions of the respec- 
tive probe transition coherences oscillating in phase with 
the incident probe field. According to this criterion, for 
A / 0, only one of the terms in Eq. ([24)) contributes to 
each probe field coherence. In contrast, for A = 0, all 
terms in Eq. ([2"4")l contribute to this response. The addi- 
tional terms lead to the cross terms discussed here. In 



principle, also other terms contribute in this case, but 
they are of higher order in the probe field Rabi frequen- 
cies and can therefore be neglected in linear response 
theory. 

How can we interpret the different contributions 
(k 6 {-1, 0, 1}) to Eq. ([M]) physically? An explicit calcu- 
lation reveals their dependence on the probe and control 
field Rabi frequencies. The control field Rabi frequen- 
cies appear for two different reasons. First, the popula- 
tions depend on the control field Rabi frequencies. But 
second, the Rabi frequencies also indicate the physical 
process described by the respective terms. The obtained 
combinations of the different Rabi frequencies lead to an 
interpretation of Rk as depicted in Figs. [5] and [6l For 
example, Rq contributing to Q21 is shown in Fig. EJa). 
This term arises from the scattering of the control field 
off of transition |1) — ► |3) and the probe field off of tran- 
sition |3) — > \2) into the probe transition \2) —> |1), which 
contributes to the magnetic response. 

We now turn to a numerical study of the cross terms. 
In Fig. we plot 6> 2 i' ^ multiplied by a factor of 7 to 

achieve unitless quantities. 0^2^ 1S n °t shown, as it is 
virtually identical for the chosen parameters. It is impor- 
tant to note that in Fig. [7J the cross terms are plotted 
over the variable a which is defined as a — — A2) = 
u>2 — h(<^3 + u>x). This new variable can be interpreted 
as the energy shift of state 2} with respect to the av- 
erage energy of |1) and |3). Hence, in a plot against a, 
effectively state \2) is moved. In this way, the probe field 
frequency remains fixed such that the multiphoton reso- 
nance condition A = 4=> u> p = lo c /2 is fulfilled for all 
values of a. For a = 0, \2) lies in the very middle of |1) 
and 1 3). 

It turns out that apart from the phases of the dipolc 
moments, the chiralities in Eqs. (Q3J) depend on the phase 
tp — 2cf) + Kf arising from the closed interaction loop, 
where K = (k c — 2k p ) is the so-called wave vector mis- 
match. For the plot in Fig. [Jj we set all involved phases 
to zero. This does not affect our final conclusion regard- 
ing the enhancement of the magnetic response, since it 
will only be based on the magnitude of the cross terms. 
If in an experiment the difference between the absolute 
field phases ip — 2(f> is not fixed, then the chiralities aver- 
age to zero. While absolute phase control is very difficult 
to achieve, relative phase control has been accomplished 
experimentally [22]. In related systems, the phase de- 
pendence of the cross terms can be made independent of 
the probe field phase, as will be discussed in Sec. El The 
observed phase dependence is a characteristic of closed 
loop systems, and has been observed in related systems 
as well [lH H(| ■ The phase can be calculated by following 
the interaction loop, and counting phases of fields that 
deexcite the atom and phases of fields that excite the 
atom throughout this loop path with opposite sign. 

Due to the dependence of the chiralities on Kf, a fur- 
ther condition for the enhancement arises, namely, that 
the so-called wave vector mismatch should vanish, i.e. 
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FIG. 5: (Color online) Interpretation of the different scattering processes that contribute to Q21 in the case of A = 0. (a) 
Rq contributes to the cross term in Q21, and thus to the chirality £he- (b) Ri contributes to the direct term, i.e. to the 
magnetic susceptibility, (c) J?_i is of higher order in either one of the probe field Rabi frequencies and therefore negligible in 
our calculation. 




FIG. 6: (Color online) Interpretation of the different scattering processes that contribute to ^32 in the case of A = 0. (a) Ro 
contributes to the direct term in £32, and thus to the electric susceptibility, (b) Ri is of higher order in either one of the probe 
field Rabi frequencies and can therefore be neglected, (c) contributes to one of the cross terms, i.e., to one of the chirality 
coefficients. 



K = 0. This condition on the relative propagation direc- 
tions of the different fields for the enhancement to take 
place can be fulfilled, for example, for co-propagating 
fields [Ij. 



C. Enhancement of the magnetic response 

We are now in the position to evaluate the magnitude 
of the magnetic response. For this, we examine the differ- 
ent contributions in Eq. (|19b[) . In particular, we compare 
the two contributions to the magnetization, which except 
for the common prefactor are given by 

Mr = d 3 2M2i £ P e l * , (38a) 
M 2 = B p g { 2 \ 0) . (38b) 

Here, Mi refers to the cross term contribution that only 
contributes at A = and M2 denotes the direct term. 
First, Figs.[2]and[7]show that in the closed loop system, 



Thus, the magnitude of the two expansion coefficients 
is comparable. The order of magnitude of the involved 
transition dipole moments can be estimated as, 

M21 ~ Ms ~ ea^ac , (40a) 
d 32 ~ ea , (40b) 

where /is is the Bohr magneton, e the elementary charge, 
ao the Bohr radius, a the fine-structure constant and c 
the vacuum speed of light. Finally, 

B p = l -Z p . (41) 

With Eqs. |@DJ) and (gTJ), we thus arrive at 

|Mi| wa-^Mal. (42) 

This means that at multiphoton resonance, the magnetic 
response of the closed-loop system is enhanced by a fac- 
tor of a -1 due to the scattering of the electric probe field 
component into the magnetic probe transition (described 
by Rq in Fig. [5]). A similar argument shows that the di- 
rect terms of the incoherently pumped and the closed 
loop system are comparable in magnitude, such that the 
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FIG. 7: (Color online) Real (red solid line) and imaginary 
part (blue dashed) of the expansion coefficient f? 2 i m the 
closed-loop system [see Eqs. flj) ]. f?^' '' is not shown, since 
it is virtually identical to the shown results. The coefficients 
correspond to cross terms and determine £he and (es, re- 
spectively. The plotted coefficients are phase-dependent; in 
this figure all phases are set to zero. The variable a = 
(A3 — A 2 )/2 = lj2 — (^3 + ui)/2 denotes the shift of the 
eigenfrequency of \2) with respect to the average frequency 
of |1) and |3). Then, for all values of a, the multiphoton 
resonance condition A = is fulfilled for a fixed probe field 
frequency uj v . 



E 





k 


-|2) 










r 



|3) 



-^-ll> 



|4) 



FIG. 8: (Color online) Example level scheme for magnetic re- 
sponse enhancement independent of the probe field frequency. 
fli and 0,2 are coherent coupling fields with frequencies u>\ and 
ui2, respectively. E and B are the electric and magnetic com- 
ponent of the probe field with frequency uj p . Here, the closed 
loop path contains an absorption and an emission of a probe 
field photon such that the multiphoton detuning A = ui\ — U2 
can be satisfied for arbitrary probe field frequencies. 



closed-loop system in multiphoton resonance also allows 
an enhancement of the magnetic response by a -1 as com- 
pared to the incoherently pumped system. 



V. DISCUSSION 

We have seen that the direct response, which describes 
the polarization [magnetization] created by the electric 



[magnetic] probe field component, is of the same order of 
magnitude in all three cases considered in Sec. IIIII Only 
at multiphoton resonance, the response of the closed- loop 
system in addition contains a term corresponding to the 
scattering of the electric probe field component into the 
magnetic probe field transition. We could show that this 
scattering effectively enhances the magnetic response by 
one inverse power of the fine structure constant a , and 
thus clearly identify the mechanism leading to the en- 
hanced magnetic response. 

Of course, the reverse process of scattering the mag- 
netic field into the electric probe field mode can also oc- 
cur. However, Eqs. (fTT)|) show that this does not lead to 
an enhancement of the electric response, since, mathe- 
matically speaking, \d\ 2 E p \ ~ a~ 1 |c?32 [132 B p \. Physi- 
cally speaking, the coupling of the magnetic probe field 
component to a magnetic transition is smaller than the 
coupling of the electric probe field component to an elec- 
tric transition by a factor of a. 

The multiphoton resonance condition restricts the 
magnetic enhancement in the system discussed here to 
a single probe field frequency. Extended closed-loop sys- 
tems can be constructed in such a way that a complete 
loop contains an excitation on the electric probe transi- 
tion and a deexcitation on the magnetic probe transition, 
or vice versa [H], EH ■ A simple example for this is shown 
in Fig. [5J In such a system, the multiphoton resonance 
condition does not depend on the frequency of the probe 
field, such that it can be fulfilled for arbitrary probe field 
frequencies by suitably choosing the coupling field fre- 
quencies. The physical interpretation identified in our 
analysis directly carries over to these extended systems. 

While the parametric enhancement by a -1 is univer- 
sal, other closed-loop systems could in principle lead to 
a modification of the relative magnitude of the differ- 
ent expansion coefficients in Eq. ([55]) . If this ratio can 
be altered favorably, then the enhancement can be even 
higher. It remains to be seen, however, whether there 
is a physical mechanism that enables one to change this 
ratio to a great extend. 

In conclusion, we have studied a mechanism for the 
enhancement of the magnetic response in atomic systems. 
We found that the studied enhancement can be traced 
back to a scattering of the electric probe field component 
into the magnetic probe field component. This scattering 
is only possible in closed-loop systems, where the applied 
laser fields form a closed interaction loop. Our analysis 
shows that the magnetic field enhancement is of order 
a -1 and occurs if the multiphoton resonance condition 
is fulfilled. 
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